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Abstract 

Models that provide experimentally testable violations of ordinary Quan- 
tum Mechanics have been recently proposed. These models are based on 
non-unitary time evolutions of density matrices that are generated by 
linear positive maps. We discuss the consequences of imposing a stronger 
condition on those maps, known as complete positivity. It turns out that 
experimental data on the neutral kaon system giving upper bounds to 
the parameters characterizing positive maps, also give bounds to those 
determining completely positive ones. 
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Since its discovery, the K-K system has been used as a testing ground for fundamental 
questions in Quantum Mechanics. At first, experiments with kaons gave precise evidence 
to violations of discrete symmetries. Nowadays, it seems possible to test the completeness 
of Quantum Mechanics itself, particularly thanks to the high luminosity of the new particle 
factories being built [1] . 

The physical intuition underlying all the recently proposed models that give rise to 
violations of Quantum Mechanics stems back to the idea that the "foam" structure of 
space-time at Planck's length produces loss of phase-coherence [2, 3]. Therefore, mixed 
states (density matrices) and no more pure states are fundamental for the description 
of physical systems, and accordingly time-evolutions transforming pure states into mixed 
ones are allowed. 

Quantum evolutions based on the standard Liouville equation 

^ = -lfi,pWl, (1) 

always preserve purity of states; indeed, if the initial density matrix p(0) is pure, i.e. 
[p(0)]2 = p(0), then for all times: 

[p{t)f = Pit) . (2) 
When the Hamiltonian H is not hermitian, Eq.(l) becomes 

dpit) 



dt 



= -iH p{t)+ip{t)HK (3) 



This happens in the phenomenological description of the decay properties of the K-K 
system (Weisskopf-Wigner approximation). We notice that Eq.(3) does not yield condition 
(2), since H ^ entails loss of probability. However, this is not loss of coherence: 
pure states remain pure, except for a (time-dependent) normalization. In fact, while (1) 
generates a time-evolution with group property, (3) produces a decay-like evolution, hence 
the group composition property holds either for positive or negative times. 

Loss of coherence shows up when further terms are added to the r.h.s. of (3). We will 
consider minimal modifications of the time evolution equation of the type 

^ = -iH Pit) + ipit) t + h[pit)] , (4) 

in such a way that the additional piece h[p] generates an evolution map rt : p(0) i-^ pit) 
that is linear and transforms density matrices into density matrices. This means that 
probability is preserved by r^: Tr(p(t)) = Tr(p(0)), and that loss of probability might 
only come from non- hermitian parts of the Hamiltonian. Further, density matrices are 
positive operators; this means that any density matrix p can be written as p = crV, for 
a suitable operator a. Asking for pit) to be again a density matrix implies that Tt has to 
be positive as a linear transformation. More generally, these minimal requests follow from 
the assumption of the existence of a scattering operator mapping linearly in-states into 
out-states [3]. 
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A stronger request to be made on the map h[p\ is that the evolution it generates, be 
not only positive, but "completely positive" [4-7]. This requirement has a precise physical 
meaning. Suppose we have two separate systems A and B, described by the density 
matrices pA and ps- If Pa evolves into a new state by means of the dynamics of A only, 
then B can not be affected by that changement, and pB remains unaltered. In other words, 
two physically uncorrelated systems can not influence each other [7]. 

More precisely, suppose the evolution in A transforms the state pA into u) \pa\ ; then, 
complete positivity amounts to requiring that the map defined by 

[pa ® Pb] =^ [pa\ ® Pb , (5) 

be positive for any separate system B. 

The simplest physical instance of a completely positive evolution is any unitary one, 
i.e. those generated by Eq.(l). More interestingly, think of the two separate systems A 
and B as forming a global system evolving under a unitary group U (t) . If p is the state of 
the compound system, the reduced states of the subsystems are given by the partial traces 
over the unwanted degrees of freedom; with obvious notations: p^ = TrB(p), ps = Tr^(p). 
Let Oa be an observable of the subsystem A. As an observable of the global system, it 
can be written as O = Oa ® 1b- Its mean value at time t follows from: 

{0{t))p = TV {pU{t)OU\t)) = Tta {Ttb [U^it)pU{t)] Oa) • (6) 

One can check that the map uJt : Pa ^ PAif) = Tfs [^^(^)P^(^)] is completely positive 
and maps density matrices into density matrices, Tryi(pyi(t)) = T^a{,Pa) [4]. 

Notice however, that the sub-dynamics ujt on A does not in general satisfy a group 
composition law Ut o ijjg — ojt+s, not even for t, s positive (semigroup composition law). 
On the other hand, an obvious physical requirement for the evolution generated by h[p\ 
in (4) is that it forms a semigroup: 

n°rs = n+s , t, s > . (7) 

The two requests, semigroup composition law and complete positivity, fixes the form of 
the generator of [5, 6], and the linear map h[p{t)] in Eq.(4) must read: 

hW)] = ( E ^]^^- ^(^) + P(^) E 4^^) + E '^At)A] , (8) 
j j j 

where the Aj are suitable operators, such that R = be bounded. 

The map h is further constrained by the request that the (von Neumann) entropy, 
S[p] = — Tr(plnp) increases: 

(9) 

dt - ^ ^ 

This corresponds to the idea that information can be lost, for example due to quantum 
gravity effects. From Eq.(8), this is most simply achieved by asking that Aj = Aj.[8\ 
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Some remarks are in order at this point. Suppose that the system B considered above 
is a heat bath. Then, there exists a canonical technique to produce a completely positive 
semigroup driving the subsystem A (a so called quantum dynamical semigroup [4-6]). It 
is based on two approximations: first one assumes a weak interaction between A and the 
reservoir B, and second that the characteristic times of B be much shorter than those of 
A. Technically, one performs a weak coupling limit (or van Hove limit) on Eq.(6) and ends 
up with a generator of the form (8) for the dynamics of A. 

In this respect, let us point out that there are two philosophically quite different 
attitudes in dealing with the evolution generated by Eq.(4) and Eq.(8): either one considers 
it as an effective description derived from a more fundamental dynamics, ore one deals with 
that evolution as fundamental in itself. An example of the latter point of view is discussed 
in [9], where the origin of the additional piece h[p] in (4) is traced back to the "foamy" 
structure of space-time, which is intrinsically unobservable, except for its effects on loss 
of quantum coherence. As a consequence, the physical time-evolution violates ordinary 
Quantum Mechanics. 

The situation of a subsystem weakly coupled to a heat bath discussed before is instead 
an example of the former point of view. In this case. Quantum Mechanics is not violated at 
the level of the global dynamics of the subsystem plus reservoir, but only by the reduced 
dynamics. Further, notice that quantum dynamical semigroups generated by (8) arise 
as effective dynamics also in situations involving environments quite different from heat 
baths: for instance, from the interaction of a microsystem with a macroscopic measuring 
apparatus [10]. Therefore, treating the modified evolution as an effective dynamics allows 
discussing in principle more general physical conditions. 

As a final remark, we stress that, in general, Eq.(4) and Eq.(8) do not respect con- 
servation laws even in presence of symmetries [9]. While from a phenomenological point 
of view this fact can be tolerated since in any case symmetries correspond to quantities 
that are conserved by the global dynamics, this is no longer true if Eq.(4) and Eq.(8) 
are regarded as fundamental. The consequences of this fact have been discussed in [11, 
12]. In any case, notice that these violations can be confined within space-time regions of 
dimensions of the order of Planck's length [13]. 

In what follows we take on the "phenomenological" attitude and consider Eq.(4) and 
Eq.(8) as effective evolution equations and apply them to study the K-K system. The 
mixing of these two particles can be modeled by means of a two-dimensional Hilbert space. 
We use the CP-eigenstates \Ki) and \K2) as an orthonormal basis: 

l^i) = 71 [l^) + 1^)] ' 1^2) = [l^) - 1^)] • (10) 

With respect to this basis a density matrix will be written as: 

Pll Pl2 



(11) 

Pl2 P22 

where * signifies complex conjugation. In the standard quantum mechanical description, 
the decay properties of the K-K system are very well captured by the Weisskopf-Wigner 
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Hamiltonian which is conventionally written as H = M—^F, where M and F are hermitian 
2x2 matrices. Wc use the convention of [14] and characterize them in terms of the complex 
parameters eg, el, appearing in the eigenstates, 

\Ks)^NsQ^^ , \Kl) = Nl(^'^^ , (12) 
and the four real parameters, m^, 75 and tulj 'Jl characterizing the eigenvalues of H: 

><s = ms- -is , \L = mL- -iL ; (13) 

Ns and Nl are normalization factors. 

Motivated by the previous considerations about complete positivity and entropy in- 
crease, we will generalize the Weisskopf-Wigner evolution equation (3) for the K-K system 
and write: 



dp{t) 
dt 



-iHp{t) + ip{t) i?t - 1 [i?p(t) + p(t) e\ +Y.^iP^t)^o ' (14) 



where the Aj are now hermitian 2x2 matrices. 

Both p and the Aj 's can be expanded in terms of Pauli matrices ai and the identity 

P = PaCra, Aj^^a'^aa, q; = 0,1,2,3. (15) 

For notational convenience, we introduce the following three vectors with components 
labeled by the index j 

ax = {a{) , a2 = (a^) , 03 = (4) . (16) 
In this way one can easily work out the action of 

h[p] = -\ \^p{t) + pit) b\^Y. ^3pii)^3 , (17) 

j 

as a 4 X 4 matrix [hap] on the column vector with components (po, Pi, P2, Ps); the entries 
hap can be expressed in terms of the scalar products of the vectors (16): 

3 

hop^hao^O, hrs ^ ■ as - 2Srs ■ ttk , r, 5 = 1,2,3. (18) 

k=l 

Notice that the Oq components of the matrices Aj do not appear and the number of free 
parameters is six: the lengths of the vectors ai, 02, 03, and the corresponding angles 
between them. Therefore, there is a minimal choice for the matrices Aj, the one for 
which: j = 1,2,3, the components Oq = 0, and ai, 02 and a^ are linearly independent 
three-dimensional real vectors. 
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Prom Eq.(18), it follows that the three-dimensional matrix [hrs] is symmetric, and 
further that its entries satisfy the following inequalities: 

hii > h22 + hss , < - {hn - ^122)^ , 

^22 > hii + h33 , < /i22 - {hii - hssY , (19) 

hss > hii + h22 , /i23 ^ ^11 - (^22 - hss)^ . 

These are direct consequence of the positivity of the norms and of the Schwartz inequalities 
for the scalar products of the three vectors ai, (12 and 03. 

Note that if one of the diagonal entries of [hrs] is zero, then the other two must be equal 
and the off-diagonal terms must vanish. For example, taking hn — 0, the four-dimensional 
matrix [/ta/j] reduces to: 



[hap] = 














\ 






















h22 







Vo 








h22 


) 



(20) 



Comparing this result with the quantum mechanical violating term h[p{t)] proposed in 
[9] , we see that there the evolution equation for the K-K system gives rise to a completely 
positive semigroup only if h[p] acts as in (20); in the notation of [9], this means that the 
three phenomenological parameters must fulfill: a = 7 and /? = 0. 

The choice made in [9] for the map h[p] is a very particular one; in fact, despite 
our requirement of complete positivity, we have three more phenomenological parameters, 
although subjected to the conditions (19). The reason for this is that in [9] it is argued 
that modifications of Quantum Mechanics coming from quantum gravity might violate 
strangeness conservation only as a second order effect. In our approach, such an attitude 
would force us to take hn = in (18), and thus to reduce to (20). 

We would like to point out that violation of the quantum mechanical time-evolution 
might be due to other effects than gravitational ones. These other effects, as briefly 
mentioned above, could be accounted for by the request of complete positivity. The larger 
number of parameters in [hap] thus available need to be fixed from the experimental 
data. We find it remarkable that the condition of complete positivity might be tested 
experimentally. 

In the remaining part of this letter we shall obtain some order of magnitude estimate 
for some of the parameters entering the matrix (18), by comparing the evolution dictated 
by Eq.(4) with experimentally measured quantities. 

Any physical property of the K-K system can be extracted from the density matrix 
p{t) by taking its trace with suitable hermitian operators. Useful observables are associated 
with the decays of the neutral kaons into 2ti or Stt states, or into semileptonic states -Kif. 

We take the approximation in which only the K meson can decay into the state tt~£'^v 
(therefore ignoring small violations of the so called AS' = AQ rule [15]). In the basis Ki, 
K2 of Eq.(lO), the operator responsible for the decay in this final state is proportional to: 

0,.= (\ \) . (21) 
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Similarly, for the decay into the state ir'^i v one has: 



1 -1 
-1 1 



(22) 



With these operators, one can form an experimentally testable observable, the so called 
CP violating charge asymmetry [9, 14]: 



m 



Tr 



P{t) (a 



Tr 



(23) 



pit) (0,+ +0,-) 

The decay of neutral kaons into pions are described by the two operators [14, 16]: 

e'2,r = ( ; n 1 , c»3. = ( : 1 . (24) 



Out of these operators one can form two new observables, the decay rates of the neutral 
kaon system into 27r and Stt: 



Tr 




Tr 


P(0)C»2. 



Tr 




Tr 


p(0)C»3,r 



(25) 



To compute the observables 6{t) and R{t), one needs to solve the evolution equations 
(4) and (8) for the density matrix (11). We shall explicitly give only the asymptotic 
expression for the solutions of those equations in the long and short time regime. A more 
detailed and complete treatment, which allows a comparison of the quantities (23) and 
(25) with the experimental data also at intermediate times, will be reported elsewhere. 

One can check that for large times, the solution of Eq.(4) with h[p\ as in (8), decays 
exponentially to 



PL = 



lerP _ has. 



h23+ihi3 
2AA 

f=* — 



fe23-»fel3 

2(AA)* 



'^L 2AA 



(26) 



while for short times, the behavior of p{t) is again exponential but with different coeffi- 
cients: 



PS I I h23-ihi3 
-r 2AA 



|e„|2 I h33_ 
\*^s\ 2Ar 



2AA 



h23+ihi3 
2(AA)* 

f + 2|^Tm 



rr. h23+ihi3 
•"•5 (AA)* 



(27) 



both these asymptotic solutions are independent from the initial conditions. In the previous 
formulas, AA — Xl — — Am + iAr/2 denotes the difference of the eigenvalues (13) of 
the Weisskopf-Wigner hamiltonian, with Am = uil — ms the mass difference of the Kl 
and Ks particles, and Ar = ^s—1l the difference of their decay widths. The two density 
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matrices pi, and ps describe mixed states, corresponding to a mixture of the conventional 
Kl and Kg neutral kaons. 

On general grounds [9, 14, 16], one expects the values of the parameters hrs to be very 
small. For example, assuming that the extra term h[p{t)\ in the evolution equation (4) is 
due to the effects of quantum gravity, a rough estimate for hrs is m\/Mp ~ 10~^° GeV, 
where Mp is the Planck mass. In writing and ps above, we have taken into account 
this fact and kept only the leading contributions in each entry of the two matrices. This 
approximation will suffice for the comparison with the experimental data to the accuracy 
we require. 

Inserting (26) and (27) into the expressions for the observables 5{t) and -R(t), one 
easily gets, to lowest order: 



and 



27^e 



— 1^ |2 



cl - 



h23 + ih 



13 



2AA 



5s = 27^e 



es + 



/i23 - ih 
2AA 



13 



(28) 



h 



33 



2Ar 
2Ar 



h23 + ih 



13 



2AA 

/i23 - ihiz 



2AA 



+ 2 



Am _ 
2—— Im 

Ar 

Am 



h 



23 



ih 



13 



Ar 



Im 



(AA)* 
h23 + ihi3 



{Axy 



(29) 



One can now express the parameters /i23, hi3 and /isa in terms of these phenomenological 
quantities. Ignoring CPT- violating effects in the Weisskopf-Wigner hamiltonian, = 
es = e, and noticing that Am ~ Ar/2, one derives from (28) and (29): 



h 



23 



h 



33 



2{ds-dL)Am, 

[4 7^e(e)-(5L + 5s)]Am, 

(i?f,-i?i) Am-2/i23 7le(e) 



(30a) 
(306) 
(30c) 



Although not all the phenomenological parameters appearing in the r.h.s. are well mea- 
sured, from the existing experimental data [17] the following order of magnitude upper 
bounds can be deduced: 



|/i23| < 10"^^GeV , l/iisl < 10"^^GeV , j/issl < 10"^° GeV . 



(31) 



To obtain an estimate on the values of the remaining three parameters in the matrix 
[hrs] one needs a more detailed treatment based on solutions of the evolution equations 
Eq.(4) and Eq.(8) valid also at intermediate times and not only in the asymptotic long- 
and short-time regions; work along these lines is in progress. Nevertheless, notice that 
inequalities (19) give some partial constraints on these parameters: \hi2\ < |/i33|/2 and 
\hii — /i22| < \h33\. Since |/i33| is very small, one expects \hi2\ to be even smaller and \hii\ 
of the same order of magnitude of |/i22|- 
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